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We present theoretical prospects for creating p-wave paired BCS states of magnetic trapped 
fermionic atoms. Based on our earlier proposal of using dc electric fields to control both the strength 
and anisotropic characteristic of atom-atom interaction and our recently completed multi-channel 
atomic collision calculations we discover that p-wave pairing with 40 K and 82 ' 84 > 86 Rb in the low 
field seeking maximum spin polarized state represent excellent choices for achieving superfluid BCS 
states; and may be realizable with current technology in laser cooling, magnetic trapping, and 
evaporative/sympathetic cooling, provided the required strong electric field can be applied. We 
also comment on the prospects of similar p-wave paired BCS states in 6 Li, and more generally on 
creating other types exotic BCS states. Our study will open a new area in the vigorous pursuit to 
create a quantum degenerate fermionic atom vapor. 

O ■ 03.75.Fi, 67.40.-w, 32.80.Pj, 42.50.Vk 

The success of atomic Bose-Einstein condensation 
(BEC) pi has induced an exponential growth of inter- 
est in the properties of ultracold dilute quantum gases. 
Of particular interest now is the physics of trapping and 
cooling of fermionic atoms M. Indeed the prospect of 
superfluidity with dilute atomic vapors has already been 
studied by several groups ||-||] . 
T , In this Letter we present a theoretical study of the 

""O ■ prospect of superfluidity in magnetic trapped fermionic 
atoms. We conduct our discussions for the three alkali 
species currently being investigated at several labs: 6 Li 
(I = 1), 40 K (/ = 4), and 82,84,86 Rb (/ = 1)2)2 ) ||_ 

This Letter is organized as follows: We briefly review the 

existing proposals for BCS states of magnetically trapped 

atoms |J|-[|, after which we present our proposal based 

on the use of external dc-electric (dc-E) fields to induce 
{vq anisotropic atom-atom interactions fa] . We then address 

\Q • its advantages and present results of our detailed multi- 
channel atomic collision calculations. We conclude with 

a discussion of other possibilities of BCS states offered 

by our proposal, notably, the analogies of the Anderson- 

Brinkman-Morel || and Balian-Werthamer |1(| states in 
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*-h p In the standard BCS theory [g_3j of superconductivity 

i-rt ' the attractive interaction between electrons of opposite 

spin and momentum causes the (Cooper) instability of 

the filled Fermi sea ground state. Physically, this effect 

>-» ' can be roughly understood as some kind of condensation 

J> , of "bosonic" Cooper pairs [|12| , |l3|| . In the weak coupling 
limit when ksT c <C Tiujd, (Ji^d is a characteristic energy 
over which an attractive interaction persists around the 

jrt ■ Fermi surface), the transition temperature T c is 



kvtTc ~ hujn exp — - - — - . (1) 

P V N(0)\V\J y ' 

N(0) is the electronic density of states at the Fermi sur- 
face, and is oc kp, the Fermi momentum. V is the two 
body scattering amplitude, assumed constant and nega- 
tive around the Fermi surface. Although derived for the 



case of a homogeneous system, as argued in several recent 
studies yH Eq. y) should also provide good estimates 
for current magnetic traps where the local density ap- 
proximation is expected to be valid. We note that the 
computation of V is a complicated and difficult many- 
body problem, and reliable values can rarely be obtained. 
For an ultracold weakly interacting dilute atom gas, the 
situation is quite different, the scattering amplitude can 
be accurately computed, and is usually dominated by the 
lowest several partial waves (at sub-mini-Kelvin temper- 
atures). The effective interaction V\ for the dominant 
1-th partial wave term is ~ 5i{kp)/kp, where 5i(kp) is 
the phase shift for the 1-th partial wave. 

In the pairing of two fermions the Pauli exclusion prin- 
ciple requires the total wavefunction to be antisymmetric. 
Therefore the orbital angular momentum I has to be odd 
for spin symmetric pairs, and even for spin asymmetric 
ones. For spherically symmetric interatomic potentials 
which vanish exponentially, 6i(k) ~ k 2l+1 at low energies 
when k — > 0. For a potential of the asymptotic form 
-1/R n (n > 4), the phase shift scales as 6i(k) ~ k 2l+1 
for I < (n - 3)/2 and as S t (k) ~ k n ~ 2 for l> (n- 3)/2. 
At a large internuclear distance R, the typical ground 
state interatomic potential is asymptotically dominated 
by the van der Waals term —Ce/R 6 , where Ce is the 
dispersion coefficient. Therefore the phase shift scales as 
Si(k) ~ k 2l+1 for / = 0, f and as S t (k) ~ k 4 for I > 2. The 
only significant term is then from the s-wave, described 
by the scattering length s sc = — lim^o 5o{k)/k. 

Based on the observation above, one realizes that 
atoms in spin symmetric states interact only through a 
vanishingly small p-wave (~ fc 3 ). For 6 Li in the triplet 
(electronic spin symmetric) state, the p-wave scattering 
length defined as a 3 c _ p = limfc_>o — 5i(k)/k 3 [Q, is about 
—35 (a.u.), corresponds to a very weak attraction. While 
the s-wave triplet scattering length a sc is enormous at 
about ~ —2160 (a.u.) due to a shallow bound state p4[ . 
In a recent paper Stoof et. al. studied the prospects for 
superfluidity of magnetically trapped 6 Li [|3|. In the max- 
imum polarized state |6), a p-wave BCS pairing occurs 
at a transition temperature of 



l c — — exp 
Kb 



2(/c_F|a sc _ p |) 3 



(2) 



too slow to be practical. For instance, with a density of 
10 12 /cm 3 , we have epjkB ~ 600 (nK) and kp\a sc ^ v \ ~ 
7x 10~ 3 @. 

Kagan et. al. suggested the trapping of atoms in mul- 
tiple hyperfine states M . They discovered that two atoms 
in the same state (|6) of 6 Li) can attract through phonon 
induced interactions (caused by density fluctuations in 
the other state, |5) of 6 Li), and pairing occurs at 



l c — -— exp 
Kb 



-13 



2kp\a s 



(3) 



where the large s-wave scattering length (for 6 Li) enters. 
Despite a dramatic enhancement (from the kp depen- 



dence), this scheme is still difficult to realize. For a 
density of about 10 12 /cm 3 per state, (a total density of 
2 x 10 12 /cm 3 ), one gets k F \a sc \ ~ 0.4. The prefactor 13 
in the exponent dictates an excessively low temperature. 
Recently Stoof et. al. B proposed an s-wave pairing 
of atoms in different hyperfine states (|5) and |6) of 6 Li) 
which occurs at 

kB \ 2kp\a s ' 

For a density of ~ 10 12 /cm 3 , T c is about a few tens of a 
nano Kelvin, potentially within experimental reach if the 
atomic decay channels are suppressed Q. There are, 
however, major obstacles towards trapping more than 
one hyperfine state simultaneously: first, the two body 
exchange loss channel is now open and the dipolar rates 
generally also increase; second, for stable BCS states, the 
pairing energies for atoms in different hyperfine states 
have to overcome the Zeeman energy deficit. An inter- 
esting proposal was also recently suggested by Modawi 
et. al. [|4| utilizing the s-wave pairing of atoms in spin 
anti-symmetric states with a result similar to Eq. (0). 

In a recent letter || we proposed to control atom-atom 
interactions at ultra-low temperatures with a dc-E. In the 
presence of a dc-E the spherical symmetry of the inter- 
acting atoms is distorted, consequently the long-range 
inter-atomic potential is modified by the addition of an 
anisotropic dipole interaction term 

V E {R) = -^P 2 {cos9), (5) 

where C E = 2£ 2 af (0)af (0), £ is the strength of the 
dc-E and a x ' (0) are the static atomic dipole polariz- 
abilities for the two atoms denoted by A and B. P 2 is 
the Legendre polynomial of order 2 and is the angle 
between the directions of the electric field and the inter- 
nuclear axis. A complete low energy scattering treatment 
reveals that the total two body (unsymmetrized) scatter- 
ing amplitude takes the form M 

/(£,£') =4tt £ t 1 ;™' {k)Y: m {k)Y Vrn ,(k'). (6) 

lm,l'm' 

The reduced T-matrix elements t\™ = lim^^o T\™ /k, 
(which has a physical meaning similar to that of the s- 
wave scattering length a sc ), are all nonvanishing asymp- 
totically a k — > 0. One can intuitively understand this 
property (T/„™ ~ k in the low energy limit) based on 
the discussion of the fc-dependence of the phase shifts of 
an asymptotic potential —1/R n . We found that the cou- 
plings between different channels are essentially given by 
V2 ~ 1/i? 3 . Thus, for each channel, the effective poten- 
tial generated by the couplings, behaves in the limit of 
large R as 1/i? 6 . Therefore the character of the asymp- 
totic i?-dependence of the total effective potential for 
each partial wave equation is decided by the diagonal 



part of the total potential. For I — the effective poten- 
tial behaves as 1/R 6 and so So ~ k. For I ^ the effective 
potentials behave as 1/R 3 and so Si ~ k. The anisotropic 
nature required new numerical techniques to ensure the 
stability of the low energy scattering calculations. De- 
tailed discussions of the complete T-matrix elements for 
all alkali metal atoms will be reported elsewhere |17J . 

Both the magnitude and sign of the t\™ are tunable 
by changing £ g] . Quite generally, we found that within 
the same I and V manifold the m = m! = is always the 
largest t\™ term, and the larger m and m! , the smaller 
t\™ . Away from resonances, smaller I and V correspond 
to larger t\™ . Such a situation allows for the general 
pairing schemes as in 3 Hc p JlQjJl^ ] . One obtain the T = 
gap equation (in the weak coupling limit), 

£/ lm I'm' k 

where the general quasi-particle excitation spectra is 



% = V4 + |C(*0I 2 , (8) 

for the homogeneous case, eg = ^M — ep is the bare 
particle energy (measured from the Fermi surface at e^). 
Such a general pairing scheme results in an anisotropic 
gap C(k) and the excitation spectra E% P, ^0|Jl2[ . In par- 
ticular when one of the t\™ is made dominant and attrac- 
tive (positive), we can simplify the analysis by assuming 

C lm (k) « A lm Y; m (k). (9) 

One then obtains the (T ^ 0) gap equation 



1 = 



47T?i 2 



^(47r)^|y ;m (fc')| 2 ^:t a nh (l(3Ez) , (10) 



k> 



with /3 = l/k B T. At T c , A im (T c ) = 0, Eq. © can be 
approximately solved to give a result similar to (fil), 

T ^£ exp (-Wki^)- (11) 

except for the substitution of V — ► — i*™. The case of a p- 
wave pairing with t\^ being the dominant (and positive) 
term is indeed true for 6 Li, 40 K, and 82 ' 84 - 86 R,b in the 
maximum polarized state. These paired superfluid states 
are the analogues of the 3 He in the A\ phase. 

Such a scheme is interesting since pairing involves only 
atoms in the same hyperfine state, yet the strength of the 
pairing can be as strongas that for a s-wave pairing in dif- 
ferent hyperfine states H||,[l8| ■ This greatly increases the 
critical temperature for the BCS pairing to occur. The 
actual values of t\™ are illustrated by our calculations 
in Figs. |, 0, and | for 6 Li, 40 K, and 8 2 >8 4 >86 Rb respec- 
tively. Although hyperfine structures are not included 



in these calculations presented. We expect they give the 
correct order of magnitude and the ^-dependence since 
collisions among atoms in the maximum polarized state 
mainly proceed along the triplet potential curves. We 
have also estimated the critical temperatures Eq. (O) for 
typical trap parameters. They are plotted in dashed lines 
with respect to the logarithmic vertical axis to the right. 
The number density at the trap center is higher than the 
10 12 /cm 3 used for previous estimates, justified since all 
inelastic processes would be slower in our scheme, the 
same reason why the BEC was originally realized in sim- 
ilar maximum polarized states. In the maximum polar- 
ized states, the spin exchange collisional loss is eliminated 
and the two-body dipolar rates are suppressed. The three 
body loss inside an electric field is a much more compli- 
cated issue currently under investigation. This rate seems 
highly reasonable since the centrifugal barrier in the p- 
wave collision channel prevents close encounters of two 
atoms at short distances. 

More generally, one may consider possibilities of pair- 
ing with t\±\ terms if resonance structures can be in- 
duced. Such resonances may also provide larger t\™ 
values at smaller dc-E fields. We are currently pursuing 
detailed atomic collision calculations including hyperfine 
structures inside a dc-E to answer these questions. In- 
side an optical trap |l£| , different Zeeman states of the 
same hyperfine spin state are degenerate, but the effect 
of the dc-E can still be induced. In such a case one 
would expect that similar to the analysis of the Anderson- 
Brinkman-Morel states and Balian-Werthamer states in 
He [9 10 1J] , one can consider all kinds of spin symme- 



3 



tries of the paired atoms |18|. This would represent a 
more interesting scenario than in the case of 3 He because 
the total spins can be much larger, e.g. for 40 K, / = 9/2 
and 7/2 respectively in the two ground hyperfine mani- 
folds jl|. ' 

As illustrated in Fig. [| for 6 Li, within the range of 
dc-E field considered, the induced p-wave interaction is 
too weak for the BCS pairing to occur at reasonable tem- 
peratures. This is a result of the small electric polariz- 
ability of 6 Li. However, for 40 K (Fig. §) and 82,84,86 Rb 
(Fig. ||), pairing can be induced to occur at temperatures 
~ hui r /kB for dc-E field of about 3,000 (kV/cm) and 
1,000 (kV/cm) respectively. This should represent a rea- 
sonable prospect to aim for within current experimental 
capabilities. The estimate of T c is based on the numbers 
for currently available cylindrical magnetic traps (with 
uj z = Xuj r , where lu z and cu r are respectively the axial and 
radial frequencies), i.e. with the width of the trap ground 
state a r = y/h/2Mu r > |ij™| ||,|§. The Fermi energy 
and the momentum are respectively €f = tuj r (6XN) 1 ' 3 
and kF = (6XN) 1 ' 6 ja r for a noninteracting fermi gas 
in these traps. The effective spatial density at the trap 
center is n(f= 0) « (AA r /6) 1/3 /(7r 2 a 3 ). Although the re- 
quired electric field is high for a dc-E configuration, such 
field strengths can be easily achieved with lasers (e.g, 
CO2 laser) in a quasi-static situation pfj . 



In conclusion, we have proposed a BCS pairing scheme 
for atoms inside an external dc-E. We have performed 
detailed atomic calculations demonstrating the strong 
anisotropic pairing interactions for polarized atoms. By 
turning the p-wave scattering amplitude to larger values, 
we can effectively tune (higher) the critical temperature 
T c . Compared with the existing proposals, our scheme 
seems to have several distinct advantages: First, the pair- 
ing only involves atoms in the maximum polarized state. 
This allows the same phase space density to be reached at 
a lower spatial density, thus lowering three-body decays. 
Second, trapping inside the maximum polarized triplet 
hyperfine states completely eliminates the large (usually 
dominant) two-body decay mechanism: those due to spin 
exchange collisions. Additionally, it suppresses the two- 
body dipolar decay processes since the final decay chan- 
nels are more restricted. Third, the applied dc-E pro- 
vides a tuning knob for optimal experimental control. 
All indications seem to be consistent with our proposal 
at the moment. We therefore encourage experimental- 
ists to attempt the implementation of the strong electric 
fields required. It may also become possible to utilize 
the low energy scattering resonances to induce dominant 
two-body attractions in d-, e-, and maybe even f-waves. 
If the fermionic gases can be confined in an optical dipole 
trap similar to those recently realized for BEC at MIT 
p9[, one may be able to create more exotic BCS states 
such as those of the A- and B- phases in 3 He, but with 
richer spin structures PJlq], 

This work is supported by the ONR grant 14-97-1-0633 
and by the NSF grant No. PHY-9722410. We want to 
thank Drs. R. Hulet and C. Green for providing triplet 
potentials for Li and Rb respectively. Most of the work 
was completed while L.Y. was a long term visitor to the 
Institute for Theoretical Physics (ITP). L.Y. would like 
to thank many enlightening discussions with participants 
of the ITP BEC workshop. In particular, he would like 
to thank Profs. J. Ho, M. Gunn, Y. Kagan, and A. J. 
Leggett for the discussion of BCS physics. He would also 
like to thank the ITP for its hospitality and the support 
of the NSF grant No. PHY94-07194 (to ITP). 



APPENDIX A: SCATTERING FROM A 
ANISOTROPIC POTENTIAL 

In this appendix we provide the essential formulation 
we developed for the scattering from the anisotropic po- 
tential containing the electric induced dipole-dipolc inter- 
action. As was shown in H, the atom-atom interaction 
potential inside a dc electric field of strength £ is 

V(R) = V (R) + Ve(R), (Al) 

where Vq is the symmetric part (for the usual collisional 
studies), and the induced anisotropic part is 

V E (R) = ~Pttcoa9), (A2) 



where Ce — 2£ 2 a^(0)af (0) is the electric induced dipole 
interaction coefficient and a x ( '(0) are the static atomic 
dipole polarizabilities of atom A and B respectively. Pi (.) 
is the Legendre polynomial of order 2 and 9 is the angle 
between the directions of the electric field and the inter- 
nuclear axis. 

With such a potential ( |Al|) , the usual partial wave ex- 
pansion has to be modified. Consistent with the standard 
definition for the scattering amplitude, we denote the in- 
cident momentum (relative motion of the two atoms) by 
k, then the scattering wave function to be computed is 

Mr) ~ e l ^ p + fe^ e ifcr , (A3) 

r 

note, here r is the relative coordinates between the two 
nucleus. The scattering is described by f(k,f) which 
(for elastic scattering) is always on the energy shell, i.e. 
the scattering momentum k' — kf. We also expand the 
incident wave according to 

e *fcV = An^Mk^Y^im^), (A4) 

lm 

with the asymptotic expansion, 

1 . 7T 

jiykr) ~ - — sin(fcr — I— -), r — * oo. (A5) 

kr 2 



Therefore 

kr 



<' A ''- ^^^sin(fer-^)Y4(fc)y im (r). (A6) 



The scattering amplitude f(k. f) can also be expanded 
onto the (complete) basis Y lm (r) to yield 

4-7T 

/ ( fc ' f ) = T^ T/m(fc)y ' m(f) - (A7) 

lm 

We then have for the scattering solution 
ri>g(f) = 0fc(r) = y £*' [ y ^ sin ( fcr " l \) + Ti m (k)e ikr - a ^] Y lm (f). (A8) 

lm 

Therefore the solution 4>t(r) can also be solved in the 
multi-channel Yi m (f) basis. We have the coupled equa- 
tion, 

Hi(f>i m (f) = J2i l '~ l (l™\V{r)\l'm'}fa, m ,, (A9) 

I'm' 

where here Hi is the effective Hamiltonian for the 1-th 
partial wave 



d?_ < 1(1 + 1) 
dr 2 



h '=-— + ± ^- ( Al °) 



Th e bo undary conditions for the scattering are as from 



Eq. (A8) 



&m ~ YC m {k) sin(fcr - I*-) + T lm (k)e ikr - il % . (All) 
Rewrite the total potential as 

V(f) = V Q + V 2 V^Y 20 (r), (A12) 

with 

w=j Jtt (A13) 

we obtain 

(lm\V(r)\l'm') = V 6 iv 5 mm , + 8 mm 'V 2 (lm\20\l'm), (A14) 

and the final multi-channel form for the scattering equa- 
tions 

WimW = Vo^m(r) + X]*''"'<HV(» : )|J / m / )^/ ro (r). (A15) 

Due to the symmetriy of the matrix element 
(Zm|20|ZW) = 0, if I + Z' + 2 =odd. We immediately 
see that the even and odd parity (1) channels decouple. 
Specifically, if we keep open only the lowest two channels, 
we obtain 

ho<Poo(r) = V <p 00 (r) - V 2 <foo(r), (A]f ,, 

hzfaoir) = (V + f VSV 2 )<p 20 (r) - V 2( f> 00 (r), { ' 



and 



Mio(r) = (V + ^V 2 )cf> w (r) ~ 3^£V 2 &o(r), 
h 3 <t> 30 (r) = (V + ^=y 2 )0 3O (r) - 3^V 2 w (r), 



MiiiO") = (V - ^gW 1± i(r) - 3^|^ 3±1 W, 
ft 3 ^3± 1 (r) = (Vo-^V 2 )03±i(r)-3yjVa0i±i(r). 

To match the boundary conditions, we further expand 
the matrix scattering element according to 

r lm (k)/k = J2 t i ;;:'(w, m ,(k). (aw) 



As discussed in detail in J8| , the nature of the asymptotic 
dipole-dipole interaction potential is such that all t\™ (k) 
become independent of k at low energies. Upon matching 
with the numerically integrated solutions, we obtain all 

The total scattering amplitde is therefore 

f(k,f) = /even ll'(k,f) + f odd u>{k,r). (A20) 

For identical particle scattering one has to worry about 
the symmetrized (anti-symmetrized) amplitude /$ (/a) 
for bosons (fermions). In the spin triplet eletronic state 
of two atoms, one has 



f s (k,r)=-j=[f(k,r) + f(k,-r)] 

= V2f cvcnW (k 1 r), (A21) 

and 

f A (k,r)=±[f(k,r)-f(k,-r)} 

= V2f oddw (k,r). (A22) 
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FIG. 1. For 6 Li with A = 1, N = 10 6 , and uo r = (2tt)400 
(Hz). We then have a r — 1.45 (/^m), e F ~ 181.7 (hui r ), 
k F « 13.5/or, and n{r = 0) = 1.36 x 10 13 /cm 3 . The solid 
lines denote t\™ in (a.u.) and refer to the left vertical scale, 
while the dashed line denotes fc_gT c (for pairing with tjjj) i n 
units of huir and refers to the (10 based) logarithmic right 
vertical scale. 



FIG. 2. The same as for Fig. 1, but now for 40 K with 

A = 0.1, N = 10 6 , and u r = (2tt)400 (Hz). We then 

have a r = 0.56 (/im), e_F ~ 84.3 (hu> r ), k F fa 9.2/a r , and 
n(f =0) =7.45 x 10 13 /cm 3 . 



FIG. 3. The same as for Fig. 2, but now for 82 - 84 . 86 Rb with 
a r ~ 0.39 (/im) and n(f = 0) ~ 2.1 x 10 14 /cm 3 . All three 
isotopes give similar results since hyperfine structures are not 
included here. The half lifetimes of the three isotopes are 1.25 
minute, 32.9, and 18.8 days respectively. 
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